Scaling Properties of ID Anderson Model with Correlated Diagonal Disorder 
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Statistical and scaling properties of the Lyapunov exponent for a tight-binding model with the 
diagonal disorder described by a dichotomic process are considered near the band edge. The effect 
of correlations on scaling properties is discussed. It is shown that correlations lead to an additional 
parameter governing the validity of single parameter scaling. 



I. INTRODUCTION 

It has been known for almost forty years that in stan- 
dard one-dimensional disordered models all states are lo- 
calized for any strength of disorder, and, that there is no 
localization transition in such systems. tlEI Formally this is 
expressed by the statement that the Lyapunov exponent 
(LE), 7, defined as 



7 = — lim — In T, 

L — >oo 2l7 



(1) 



where T is the transmission coefficient ^through the sys- 
tem of length L, is always positivefl While it may 
seem that there is nothing more to be said about one- 
dimensional models, their localization properties have 
recently attracted a great deal of attention. This re- 
newed interest has been concentrated in two areas. The 
first area includes studies of unusual models that demon- 
strate the presence of extended states. It was noted first 
in Ref. || that introducing correlations in the statistical 
properties of random site energies in the Anderson model, 
one can create extended states in the one-dimensional 
Anderson model at certain values of energy. Later on, 
the authors of Ref. ^ showed that regions of extended 
states exist in ID systems with long range correlations 
of the random potential. The authors of Ref. |^ demon- 
strated that vanishing LE can be obtained in an arbi- 
trarily chosen spectral region by selecting a special form 
of the correlation function of the random potential. It 
should be noted, however, that this result was obtained 
in Ref. ^ only in the second order of the weak disorder 
expansion. Taking into account the next terms of the 
expansion makes LE positive again. Q 

The second area of research interest in the field of one- 
dimensional localization is focused upon scaling and sta- 
tistical properties of finite size systems. The single pa- 
rameter scaling put forward in Ref. |l| is a cornerstone of 
the current approach to the Anderson localization, but its 
interpretation in the presence of non-self-averaging fluc- 
tuations of conductivity was the subject of long debate. 
Eventually it was understood that the SPS hypothesis 
in Anderson localization means that the entire distribu- 
tion function of conductivity (or transmission), and not 
individual moments, must be parameterized by a single 
parameter.l3i3 This distribution iii one-dimensional sys- 
tems has been intensively studiedB, and it was established 
that the bulk of the distribution for sufhciently long sys- 
tems has a log-normal form. It was also understood start- 



ing with Ref. 1^ that under certain circumstances the two 
parameters of the distribution, the average value of the 
finite size LE, 7 = — lnT/(2L), which coincides with its 
limiting value, 7 — (7), and the variance, = var(f)), 
are related to each other in a universal way 



r- — - 1 
7 



(2) 



This relation reduces two parameters of the distribution 
to only one, and provides, therefore, a justification and 
interpretation for SPS. Conditions, under which Eq. (H) 
holds, however, have been established only recently, when 
the authors of Ref. |l^ showed that the validity of SPS 
is controlled by a new macroscopic length. Is, defined in 
terms of the integral density of states, N{E): 



l-^ = sin(7r7V(£;)). 



(3) 



where E is energy. It was shown that SPS holds when 
K — (7/3)^^ > I7 and fails in the opposite case. Recently, 
the authors of the present paper showed that parameter 
K plays a more important role than just establishing the 
criterion for SPS. It was found that in the non-SPS region 
the function r defined in Eq. and the similar dimen- 
sionless combination of the third moment with the length 
L and the Lyapunov exponent 7 are functions of this sin- 
gle parameter .EJ In the recent paper Ref. [ij, higher mo- 
ments of the distribution of the Lyapunov exponent were 
calculated analytically for a quantum particle in a poten- 
tial described by the Gaussian delta-correlated random 
function. It was found that all moments of the distribu- 
tion for this model can be expressed as functions of the 
localization length and the parameter E/D'^/'^ , where 
D characterizes the strength of the potential. For the 
Gaussian white noispjDotential the parameter k is also a 
function of thus in the particular case of the 

Gaussian white noise potential these two parameters are 
equivalent. 

In the present paper, we combine the two discussed 
areas and focus upon correlation induced changes in the 
probability distribution function of the Lyapunov expo- 
nent. In order to study this problem, we consider a tight- 
binding model with the random potential described by a 
zero-mean dichotomic process. This process is character- 
ized by an exponential correlation function with a cor- 
relation radius, which is a free parameter of the model. 
There is no localization-delocalization transition in this 
model, and we are interested in effects of the correla- 
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tion radius on scaling and statistical properties of con- 
ductance in the regime of strongly locahzed states. This 
model allows for an approximate analytical treatment in 
the spectral region, where the correlation radius becomes 
a dominant length. Combining analytical and numerical 
calculations we show that the presence of this additional 
length significantly changes the transition between SPS 
and non-SPS regions, and that in the non-SPS region it 
results in a new scaling behavior. 



II. QUASICLASSICAL APPROXIMATION FOR 
LYAPUNOV EXPONENT 



We consider the tight binding model, described by the 
equation of motion 



V'n+l + + {VCn ~ E)xjj,, = 0, 



(4) 



where V describes the strength of the random potential. 
The potential randomly takes one of two values ztV de- 
pending upon whether the dichotomic random variable 
is equal to 1 or —1. The assumed value remains con- 
stant within regions of random lengths (segments), which 
are distributed according to the discrete analog of the 
Poisson distribution. The average length of such regions 
is 



L = coth 



1 



where Tc is the correlation radius defined by 



(5) 



(CnC™> = e-l"-™l/'■^ (6) 

The finite size LE is defined through the norm of the 
transfer matrix: 



^{E) - ilog||rjv...Ti| 



(7) 



where T„ are transfer matrices which, in our case, are 
equal to either T+, or T_ depending upon the value of 
Cn, where T± are 



T± = 



E±V -1 



1 







(8) 



and —2 — V < E < —2 + V, extended states exist in a ho- 
mogeneous system with the site potentials equal to —V 
or +V, respectively. Finally the third region, 2 + V < \E\ 
corresponds to the energies, where extended states would 
not arise in homogeneous systems with either of the two 
possible values of the potential. Obviously, no states can 
arise in this region also in the system with the random 
distribution of the potential. 

While all states in this model are localized, different 
spectral regions still have different transport properties. 
Transport in the first region is characterized by multiple 
scattering of propagating modes, and the localization in 
this region results from the interference of the scattered 
waves. Transport in the second region can be described 
as tunneling through states, arising within the segments 
of the system, where the potential produces potential 
wells between barriers formed by segments with the op- 
posite value of the potential. If the average length of the 
"barrier" regions is larger than the penetration length 
under the barrier, the transport is mostly determined by 
the under-barrier tunneling, and the interference effects 
related to phase changes of the wave functions inside po- 
tential wells become less important. We are mostly inter- 
ested in this region, where correlations are expected to 
result in the most non-trivial effects. One could expect 
that the transport in this region can be described within 
the quasi-classical approximation, and we will show be- 
low that this is indeed the case. 

In the context of our problem, this approximation 
means neglecting commutators between transfer matri- 
ces at different sites. Introducing LE in a homogeneous 
system with the potential U according to 



70 (S; C/) =Re [Ao(^;;7)], 



(9) 



where 



Ao = cosh 



E-U 



(10) 



we can write a quasiclassical expression for LE for Eq. 
in the following form 



1 ^ 



(11) 



To consider the spectrum of this system let us note 
that an addition of a constant potential ±t7 shifts a con- 
duction band, —2 < i5 < 2, of the homogeneous tight- 
binding system with no on-site potential by the value of 
the potential: —2 + U < E < 2 + U . Correspondingly, 
the spectrum of the random system under consideration 
can be divided into three regions with the qualitatively 
different spectral and transport properties. One region, 
< \E\ < 2 — V, corresponds to the common part of 
the conduction bands of homogeneous systems with the 
site potentials equal to either +V or —V. In the second 
spectral region, consisting of energies 2 — V<E<2 + V 



Using the fact that the Cn = ±1, Eq- ( pT| ) can be rewrit- 
ten in the form 



N 



1 = la + J^Cn 



(12) 



n=l 



where 



la = -ME;V)+ME;-V)], 
6 = ^[j,{E;V)-ME--V)]. (13) 
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The averaged finite-length LE, 7 = (7), is determined 
by 7a, while its variance cr^ depends upon S: 

= b') - {if = ^^Ic 

In the first spectral region this approximation yields van- 
ishing LE, which is quite understandable, because inter- 
ference effects responsible for the localization in this re- 
gion are completely neglected. 

Eq. (|l^) can be formally obtained in the following way. 
The relation between LE and the eigenvalues of the trans- 
fer matrix allows us to write 



2cosh(7L) =Tr {Tn,i) ■ 



(14) 



where the transfer matrix T/y 1 consists of sequences of 



Tn.i = Tn . . .Ti 



rpns rjin2 rpTli 



(15) 



Powers Uk of transfer matrices T± represent the lengths 
of the regions where the potential remains constant. The 
matrices T± can be diagonalized = R±S!IR^^ , where 



SI 



e-^°(±^) 



(16) 



Thus, introducing ai2 = R^^R- and a2i — RZ^R+, 
and noting that for long enough samples the effect of 
matrices R± appearing at the end of the structure on LE 
is negligible, we obtain that LE can be found from Eq. 
(O) with the transfer matrix replaced by 



T, 



N.l 



■ a2io.ai2S_. 



(17) 



The last step in the procedure is to commute, for ex- 
ample, S+ and ai2, with the purpose of collecting all 
S± together. If one begins this rearrangement from the 
right side of Eq. (|l^, each of such permutations pro- 
duces two terms. One of them consists of the product 
of the partially ordered and the remaining non-ordered 
parts, while the second one contains commutators of 5+ 
and ai2. Repeating the procedure leads to the recurrent 
equation for the part of the transfer matrix containing 
only products of commuting diagonal matrices S± 



T-i 



N.l 



o 



N.l 



E 



(18) 



where we have used the fact that 012021 = 1, and intro- 
duced the commutator r„ = [5*", 012]. Notation C'pj^_ J ^1 
stays for the product of all S± between sites 1, and Pk-i- 
On,i = ■ ■ ■ S^^^S^^ , where pk = V^ f-i ni. Retaining the 
first term only and using Eq. (|14D we arrive at Eq. (pi]). 

Eq. (18) can be used to obtain a representation for the 
transfer matrix in terms of the polynomial in powers of 
the parameter 6, which enters the expression for commu- 
tators r„, and can be considered as a small parameter 



in the case of weak disorder. However, the resulting ex- 
pression for the transfer matrix cannot be interpreted as 
a weak disorder expansion for LE. The reason for this is 
the fact that the power of S is determined by the number 
of segments with different values of the potential, and the 
whole polynomial is actually an expansion in terms of the 
number of "jumps" between V and —V. For instance, the 
term with results from configurations with only one 
such jump, and there are only N such configurations. In 
the limit N 00 they give zero correction to the LE (in 
addition to the exponentially small probability of such 
configurations). Thus, we can see that for long enough 
samples only terms of the order of A/" ~ N/2 make a 
significant contribution to the LE because the number of 
such terms is proportional to the binomial coefficient Cjy 
and, therefore they give a linear in N contribution to the 
log||r||. The actual number of configurations contribut- 
ing to LE depends upon the average length of the seg- 
ments with constant potential; and the number of jumps, 
corresponding to the optimal configurations can be esti- 
mated as N/lc- It is clear, therefore, that corrections to 
the first term of Eq. ( p!8| ) decrease with increasing l^- We 
thus expect that the approximation should work when 
hoc/lc <C 1. It lacjttially becomes exact in the limit of in- 
finite Tc and y.t3llil This can be understood from the fol- 
lowing arguments. Deviations from Eq. ([l^ ) are caused 
by the interference of the scattered waves. But in the sec- 
ond spectral region such scattered waves have to tunnel 
through the barriers, where the magnitude of the scat- 
tered waves decreases and the interference is destroyed. 
Obviously, the effectiveness of the destruction depends 
on the width of the potential barriers, and the rate of de- 
cay of the wave functions in them. These characteristics 
are determined by Ic and the magnitude of the potential 
respectively. 

In the spectral region, 2 — V < \E\ < 2 + V only 
7o(i?; V) in Eq. ( p^ is different from zero, and we have 
for cr^ 



(19) 



This expression is obviously different from SPS Eq. (||), 
and is valid when jlc 3> 1. In this case, scaling described 
by the function t(k) is no longer valid, but Eq. ( [is| ) sug- 
gests that a new type of scaling appears. It can be conve- 
niently described with the help of a new scaling function. 



i'lc 



(20) 



The transition between the two scaling regimes is con- 
trolled by the parameter ^Ic- If the correlation radius 
is so small that this parameter remains much less than 
unity for the entire second spectral region, the regular 
scaling behavior expressed by t(k) persists. The behav- 
ior expressed by Eq. ( pO| ) occurs when 7/c is greater than 
unity for the entire region 2 — V < \E\ < 2 + V. More de- 
tailed analysis of the transition between the two types of 
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Energy Energy 

FIG. 1: Dependence of the mean of the Lyapunov exponent 
on energy: numerical result (dotted line) and theory (sohd 
line), for potential V — 0.4 and correlation radius rc — 10. 
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FIG. 2: t(k) for different distributions of the potential: the 
uniform distribution (empty circles) and the dichotomic pro- 
cess (filled circles) when is the smallest length of the system. 



behavior requires a numerical approach, which we discuss 
in the next section of the paper. 

In the third spectral region, \E\ > 2 + V, S decreases 
with energy. This behavior can be approximately de- 
scribed as 

S^^ ^ (21) 

(V2/H~T+i) 

for < |e| - K 1, where e = (2 - \E\)/U and 

^'^(l) (log|C/|+log|e|)2 (2^) 

for |e| > 1. 

III. SCALING PROPERTIES OF LE: 
NUMERICAL SIMULATIONS 

For the numerical analysis, we calculate the LE itera- 
tively using Eq. (|^) and the usual technique of renormaL 
ization of the resultant vector after every 10 iterations.Ej 
The length Ig was calculated according to Eq. (^, where 
the integrated density of states was obtained using the 
phase formalisma with consequent averaging over all re- 
alizations. To investigate scaling properties we kept the 
length of chains much larger than the correlation length, 
Ic, and the localization length, lioc, for all magnitudes 
of the random potential U and values of the energy E 
in the vicinity of the second spectral region. Statistics 
were collected from 40000 realizations. Figs, la and b 
present results of numerical computations of the Lya- 
punov exponent. The numerical value of the correlation 
radius in these figures is such that Ic > hoc for the en- 
tire 2 ~ V < \E\ < 2 + V spectral region. A comparison 



of numerical results and the approximation ( [III ) is pro- 
vided for both average LE and its variance. For energies 
inside the first region this approximation gives zero, as 
discussed above. In the second region, one can see that 
Eq. (|l^) provides a reasonably good approximation for 7, 
and even better agreement for a^. In the third spectral 
region the agreement between Eq. (|l^) and the numerical 
calculations becomes perfect. 

We also studied numerically the scaling properties of 
our model. Let us recall that for systems without cor- 
relations SPS is controlled by the ratio k = hoc/ Is of 
two macroscopic lengths existing in the system. In the 
presence of the correlated disorder the situation is more 
complicated because a new length related to the correla- 
tion radius should be taken into account. In our model 
this length is Ic— the average length of regions with con- 
stant potential. If Ic is the smallest length in the system, 
Ic < {hoc, Is), the presence of correlations does not af- 
fect the scaling properties. Computer simulations show 
that the scaling parameter r behaves in this case qual- 
itatively similar to a system with uniformly distributed 
non-correlated disorder. In Fig. ||the dependence of r on 
K is shown for both a tight-binding model with uniformly 
distributed {—U < Vn < U) on-site potential (with U 
ranging from 0.09 to 0.145) and the model under consid- 
eration. For the latter, potentials {V = 0.001 - 0.004) 
and correlation radii {rc — 2 and 3) where chosen such as 
to have Ic the smallest length up to the genuine spectral 
boundary. 

However, when Ic > hoc, the behavior of changes 
significantly, and the character of this change depends 
upon the position of the spectral point, Ec, at which 
Ic — hoc{Ec)- If Ic is so large that Ec falls in the first 
spectral region \E\ < 2—V, the scaling function t remains 
essentially equal to unity up to the very small vicinity 
of the boundaries of this interval. Inside a very small 
neighborhood of points \E\ = 2 — V any scaling behavior 
disappears, but immediately outside of this neighborhood 
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FIG. 3: The crossover between two scaling regimes, r ~ 1 
(filled circles) and 5 ~ 1 (empty circles), as energy variable 
e = {E + 2) /V changes, is shown for a collection of different 
values of parameters {V = 0.2 — 1.2 and = 5 — 250) meeting 
the requirement 0^1. 

inside the second spectral interval 2 — V < \E\ < 2 + 
V a new scaling presented by Eq. ( pO| ) emerges. This 
situation is presented in Fig. |3[ In the third spectral 
region the function S decreases in agreement with the 
analytical formulas Eq. and Eq. (|2|) . It can be seen 
that the scaling in terms of variable e = (2 — \E\)/U 
predicted by Eq. ( ]2l] ) persists over quite a wide interval 
of energies. 

In the case when the energy Ec belongs to the interior 
of the second spectral region, the scaling properties are 
determined by the parameter k. If K(i?c) ^ 1 we re- 
turn to the situation of the non-correlated disorder and 
if k{Ec) > 1 the system does not exhibit any scaling be- 
havior in the second spectral region. The destruction of 
scaling in this case may also be illustrated by the behavior 
of t{k) for different values of Ic- When Ic increases from 
the values corresponding to the non-correlated limit, the 



function t{k) starts changing as shown in the insert in 
Fig. 2 where the dependence of r is depicted for three 
different sets of parameters of the random potential U 
and re- It deviates from the universal (SPS) dependence 
at larger values of n with the increasing magnitude of 
the deviation depending upon both Vc and U. It demon- 
strates in that way the absence of scaling properties. 



IV. CONCLUSION 

We considered statistical characteristics and scal- 
ing properties of the Lyapunov exponent in a random 
medium with correlated disorder, using a tight-binding 
model with a random potential described by a dichotomic 
process. A simple theoretical approach for estimation of 
the statistics of the Lyapunov exponent was suggested. 
Qualitative analysis based on this model allowed us to 
estimate the effect of correlations on scaling properties, 
and express this effect in terms of the ratio of the lo- 
calization length and the correlation length, 9 = Ic/hoc- 
When the correlation radius is much smaller than the lo- 
calization length {6 1), the system behaves as though 
correlations are absent and demonstrates single param- 
eter scaling, a^Llioc ^ 1, whose validity is governed by 
the ratio hoc/ Is- In the opposite case, when 0^1, the 
crossover to another scaling, a^Llf^^/lc ~ 1, specific for 
the model under consideration, takes place at values of 
energy close to the band gap. 
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